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ABSTRACT 

We calculate the number of damped Lya absorbers expected in various 
popular cosmological models as a function of redshift and compare our 
predictions with observed abundances. The Press-Schechter formalism is used 
to obtain the distribution of halos with circular velocity in different cosmologies, 
and we calibrate the relation between circular velocity and absorption cross- 
section using detailed gas dynamical simulations of a "standard" cold dark 
matter (CDM) model. Because of this calibration, our approach makes more 
reahstic assumptions about the absorption properties of collapsed objects 
than previous, analytic calculations of the damped Lya abundance. CDM 
models with Qq = 1, Hq = 50 kms~^ Mpc~^, baryon density = 0.05, and 
scale-invariant primeval fluctuations reproduce the observed incidence and 
redshift evolution of damped Lya absorption to within observational uncertainty, 
for both COBE normalization {as = 1.2) and a lower normalization (ag = 0.7) 
that better matches the observed cluster abundance at z = 0. A tilted (n = 0.8, 
cTg = 0.7) CDM model tends to underproduce absorption, especially at z = 4. 
With COBE normalization, a CDM model with Qq = 0.4, = 0.6 gives an 
acceptable fit to the observed absorption; an open CDM model is marginally 
acceptable if Qq > 0.4 and strongly inconsistent with the z = 4 data if = 0.3. 
Mixed dark matter models tend not to produce sufficient absorption, being 
roughly comparable to tilted CDM models if Q^, = 0.2 and failing drastically if 

= 0.3. 
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1. Introduction 

Damped Lya (DLA) absorbers — hydrogen absorption systems with neutral column 
density A^hi > 10^°'^ cm~^ — are important tracers of high density structure in the early 
universe. They are far more common than quasars or luminous radio galaxies, and various 
lines of evidence, both coincidences of numbers and properties of observed systems, indicate 
that DLA systems arc the high-rcdshift progenitors of "normal" present-day galaxies (see, 
e.g., Wolfe et al. 1995; Djorgovski et al. 1996; Fontana et al. 1996; Prochaska & Wolfe 1997). 
In this Letter, we use a semi-analytic method, calibrated against numerical simulations, 
to test the ability of current theories of structure formation to reproduce the observed 
abundance of DLA absorbers at 2 < 2; < 4. 

The first attempts to employ DLA absorption as a test of cosmological scenarios used 
the Press- Schechter (1974) formalism to predict the abundance of dark matter halos and 
adopted simplifying, usually conservative assumptions about the amount of gas within 
these halos that would cool and become neutral (Kauffmann & Chariot 1994; Mo Sz 
Miralda-Escude 1994; for recent analyses along similar lines see Liddle et al. 1996ab). 
Ma & Bertschinger (1994) and Klypin et al. (1995) proceeded in a similar fashion, but 
computed the halo mass function from dissipationless N-body simulations. Kauffmann 
(1996) presented a more sophisticated semi-analytic treatment of DLA systems and disk 
galaxy formation in the standard cold dark matter (CDM) model, incorporating more 
realistic assumptions about the cooling and settling of gas within halos and subsequent star 
formation. 

A fundamentally different approach was taken by Katz ct al. (1996; hereafter KWHM), 
who used numerical simulations incorporating dark matter and gas in the presence of a 
photoionizing background to show that damped Lya systems originate naturally from initial 
conditions similar to those which are believed to have produced the observed large-scale 
structure. In the KWHM simulations, damped absorption results from concentrations of 
dense gas embedded in more extended and more massive dark matter halos. The masses 
and abundances of these objects are comparable to those of galaxies, supporting the view 
that damped Lya systems are a natural by-product of galaxy formation. Lower column 
density, Lya forest absorption arises in more diffuse, highly photoionized structures, and 
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the combined results of KWHM and Hernquist et al. (1996) show that a simulation of the 
standard CDM model reproduces the observed distribution of neutral hydrogen column 
densities quite well over most of the range 10^'^ cm~^ < A^hi < 10^^ cm~^. 

Gardner et al. (1997; hereafter GKHW) improved on KWHM by using a semi-analytic 
method to correct for damped absorption in halos below the KWHM resolution limit. Using 
KWHM's simulation and higher resolution simulations focused on individual, collapsing 
halos, GKHW compute the mean projected area a{vc,z) over which a halo of circular 
velocity Vc produces damped Lya absorption at redshift z. They convolve this function with 
the Press-Schechter formula for the abundance of halos as a function of Vc to obtain the total 
incidence of damped absorption in all halos. The correction for unresolved objects raises 
the KWHM predictions by roughly a factor of two, bringing them into good agreement 
with the observed abundance of damped absorbers, though the predicted number of Lyman 
hmit systems remains low by about a factor of three. GKHW also show that star formation 
(with the algorithm of Katz, Weinberg, & Hernquist 1996) has little effect on the predicted 
number of damped systems. 

The method adopted in KWHM and GKHW can be apphed to other cosmological 
scenarios, and we believe that numerical simulations with gas dynamics will eventually 
allow the most robust confrontations between predicted and observed properties of DLA 
absorbers. However, the fully numerical approach is computationally very expensive, and 
it will be some time before it can be applied to a broad parameter space of cosmological 
models. This seems like an opportune moment to revisit the semi-analytic approach 
for several reasons. First, the theoretical ground has shifted somewhat since the first 
investigations, prompted in part by improved analyses of the COBE data (and the 4-year 
COBE data themselves) and in part by further studies of the theoretical models, especially 
the emergence of open-bubble inflation models as a viable and predictive setting for 
structure formation (e.g., Ratra & Peebles 1994; Bucher, Goldhaber, & Turok 1995; 
Yamamoto, Sasaki, & Tanaka 1995). Second, new data have improved the observational 
constraints, especially at z ^ A (Wolfe et al. 1995; Storrie-Lombardi et al. 1996). Third, 
and most significant from the point of view of this paper, the GKHW results provide a way 
to calibrate the most uncertain element of the semi-analytic method, the relation between 
halo circular velocity and the cross-section for damped absorption. It seems plausible that 
this relation is driven primarily by "local" physics within the halo, and for this paper we 
will therefore assume that the function a{vc, z) found by GKHW for the standard GDM 
model also applies to other cosmological scenarios. With this assumption — uncertain, 
but probably the best approximation available until alternative models arc investigated 
numerically — we can use the GKHW method to compute the abundance of DLA absorbers 
in any specified cosmology with Gaussian initial conditions. We describe this method in 
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more detail in §2, then present our results and discuss their implications in §3. 

2. Method 

KWHM used TreeSPH (Hernquist & Katz 1989) to evolve to redshift z = 2 a periodic 
cube 22.22 comoving Mpc on a side containing dark matter and gas. The initial conditions 
were drawn from a CDM power spectrum with ag = 0.7, Qq = 1, h = 0.5[], and = 0.05, 
where ag is the present-day rms mass fluctuation within a sphere of radius 8 h^-^ Mpc, 
and Qf) and Qq are the present-day fractions of closure density in the form of baryons and 
matter respectively. The simulation included a spatially uniform background radiation field 
of intensity J(z/) = Jq{vq/ v)F{z)^ where z/q is the Lyman-limit frequency, Jq = 10~^^ erg 
s~^ cm-2 sr-i Hz-\ and F{z) = if z > 6, F{z) = 4/(1 + z) if 3 < z < Q, and F{z) = 1 
if 2 < z < 3. More detailed descriptions of the simulation and methods can be found in 
GKHW, KWHM, and Katz et al. (1996). 

Given their limited particle number, KWHM were able to resolve halos only down 
to a mass corresponding to a virial circular velocity of Vc ~ 100 km s~^ (about IO^^Mq 
at z = 2; see, also Weinberg, Hernquist, & Katz 1997). Since Quinn, Katz, & Efstathiou 
(1996) find that halos with Vc > 37 km s~^ are capable of hosting damped Lya absorbers, 
it is clear that KWHM's simulation systematically underestimates the true prediction of 
this CDM model for the incidence of damped absorption. GKHW combine the background 
radiation field of KWHM with the initial conditions of Quinn et al. (1996) to simulate 
small regions with sufficient detail to resolve halos with circular velocities below 50 km 
s~^. The high resolution volumes are too small to calculate the number of damped systems 
(they were chosen to lie in regions where small halos would form; see Quinn et al. 1996), 
but by using the halos from both the high resolution and KWHM simulations, GKHW 
estimate a{vc,z), the cross-section for a halo of circular velocity Vc to produce damped 
Lya absorption at redshift z. GKHW show that a{vc,z) can be fitted by a power law in 
Vc, with an exponential cutoff to reproduce the absence of absorbers in halos with Vc < 37 
km s~^ (Quinn et al. 1996). The fitting procedure is described in GKHW, which also lists 
parameters of the fitted relation at 2; = 2, 3, and 4. 

According to Press-Schechter theory, the number density N{M, z) of collapsed objects 
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in the mass range M — ^ M + dM at redshift z is 

where po is the mean comoving mass density and -R/ is the Gaussian filter radius given 
by M = (27r)^/^po-R/- The parameters (Tq, 7 and -R,f are related to moments of the power 
spectrum and are defined in Bardeen et al. (1986). As detailed in GKHW, multiplying 
N{M,z) by a{vc, z) and integrating from Vc to infinity yields the number of damped 
absorbers per unit redshift residing in halos of circular velocity greater than Vc- 

(if roo 

n(z,Vc) = — N{M',z)a{v'^,z)dM', (2) 

dz Jm{vc) 

where v'^ is the circular velocity at the virial radius of a halo of virial mass M', and r is 
comoving distance. With Vc = 37 km s~^, the lower limit for damped absorption found by 
Quinn et al. (1996), equation (|^) yields the total incidence of DLA absorption n{z), defined 
to be the number of absorption systems per unit redshift with A'^hi > 10^°'^ cm~^. In this 
paper, we adopt the a{vc,z) relations found by GKHW for the standard CDM model but 
combine them with the halo abundance N{M, z) predicted for other cosmological scenarios. 

For all the Press- Schechter predictions we use a Gaussian filter with Sc = 1.69. This 
6c is higher than that conventionally used with a Gaussian filter because of the way 
that we identify halos and assign them circular velocities (see GKHW), which we believe 
most closely corresponds to the conception of halos in the Press-Schechter formalism. As 
explained in GKHW, this value of 6c fits the numerically derived halo mass function of 
KWHM quite well. Our full procedure correctly predicts n{z,Vc) in the KWHM simulations 
for damped systems in halos above the Vc = 100 km s~^ resolution limit. In our calculations 
of N{M,z), we use the transfer functions of Bardeen et al. (1986) for cold dark matter 
models and Ma (1996) for mixed dark matter models. 

3. Results 

Table [l| compares the predicted numbers of DLA absorbers to the observational results 
reported by Storrie-Lombardi et al. (1996). These results are illustrated in Figure |l], which 
also includes the observational data of Wolfe et al. (1995). To aid the interpretation of this 
Figure, we show in Figure |^ the comoving number density of halos with Vc > 37 km s~^ 
obtained by integrating equation ([l|) for each model. Starting at high redshift, the halo 
abundance first rises as rms fluctuations on the Vc ~ 37 km s~^ mass scale go nonlinear, 
then declines as the first generation of halos above this threshold merges into a smaller 
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number of more massive halos. For some of our models, this decline has already begun by 
z = 5. 

The model producing the most DLA absorption and the best agreement with the z = 4 
data is CCDM: COBE-normalized, f2o = 1 CDM with scale-invariant {n = 1) primeval 
fluctuations and h = 0.5. However, this model is known to produce excessively massive 
clusters at z = 0. SCDM, the "standard" CDM model used in GKHW, has an identical 
power spectrum shape but is normalized to produce more reasonable (still somewhat high) 
cluster masses, while disagreeing with the COBE results. In this case, n{z) is lower due 
to the decrease in power (c/. Figure but the model still matches the DLA observations 
within current uncertainties. It is possible to construct an f2o = 1 CDM model that 
simultaneously fits cluster masses and COBE data by further adjustments, such as tilting 
the primeval power spectrum to n < 1 and raising the baryon fraction. White et al. 
(1996) give an example of such a model, which we designate TCDM: with n = 0.8 and 
Qb = 0.1, COBE normalization implies ag = 0.7 just as for our SCDM model. However, 
the tilt and high Qh reduce the amount of small-scale power relative to SCDM, reducing 
the halo abundance at 2; > 2 (Figure ^ and producing an n{z) significantly below the 
Storrie-Lombardi et al. (1996) data for z = 3 and z = 4. There are other combinations of 
n, Qb, and h that can match COBE and still have CDM, Qq = 1, and ag ~ 0.7, but since 
the necessary changes all reduce the amount of power on small scales, the problem shown 
here is likely to be generic to all such models. 

It is important to note that TCDM is the only model we analyze that has flbh'^ 
different from 0.0125, and we have not allowed the increased baryon density to change the 
DLA cross-section in our treatment here. It is not difficult to imagine that the increase 
in baryon density could alter the absorption produced in high-redshift halos. Larger Qb 
implies an increased gas supply for a halo of a given mass, which could simply lead to a 
higher incidence of damped absorption. On the other hand, it could also lead to more 
rapid cooling and more dense concentrations of cold collapsed gas, thereby decreasing n{z). 
Consequently, our results for this model are not robust enough to rule it out until the 
influence of Qb on n{z) has been more thoroughly investigated. 

CDM models with a cosmological constant retain the flat universe preferred by 
inflationary models while decreasing the density of matter, permitting more reasonable 
cluster masses with a COBE normalization. They also match the shape of the power 
spectrum on galaxy cluster scales more closely than SCDM, and they avoid the severe 
problems of JIq = 1 models concerning the age of the universe and the cluster baryon 
fraction. Our LCDM model, COBE-normalized with Qq = 0.4 and h = 0.65, generally falls 
within the la errors of the Storrie-Lombardi data except ai z = 4. Many of the systems 
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Fig. 1. — The predicted number n{z) of DLA systems per unit redshift for various 
cosmological models vs. observational data. The solid error crosses (Icr) are reproduced 
from Storrie-Lombardi et al. (1996), with the vertical dashed error bars indicating 2a error. 
The dotted error crosses (la) are taken from Wolfe et al. (1995). The heavy curves connecting 
2 = 2,3,4 are values calculated for each model. The model names correspond to the names 
in Table |l], where their parameters are listed. 
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0-8 
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Other 


z = 2 


z = 3 


Z = 4: 


CCDM'^ 


1.2 


1 





0.5 




0.203 


0.310 


0.397 


SCDM 


0.7 


1 





0.5 




0.179 


0.174 


0.194 


TCDM^ 


0.7 


1 





0.5 


n = 0.8 


0.103 


0.0811 


0.0530 


LCDM^ 


1.1 


0.4 


0.6 


0.65 




0.119 


0.146 


0.155 


LCDM2'= 


0.79 


0.4 


0.6 


0.65 




0.0926 


0.0889 


0.0730 


OCDM^ 


0.65 


0.4 





0.65 




0.0916 


0.0926 


0.0836 


0CDM2'^ 


0.45 


0.3 





0.65 




0.0293 


0.0206 


0.0128 


MDM^ 


0.782 


1 





0.5 


= 0.3 


0.0246 


0.00692 


0.00243 


MDM2^ 


0.808 


1 





0.5 


= 0.2 


0.0875 


0.0545 


0.0432 










1.5 < 


z < 2.0 


2.0 < 2 < 3.0 


3.0 < ^ < 3.5 


3.5 < ;z < 4.7 


Observed ^ 








0.14 ±0.073 


0.18 ±0.039 


0.21 ±0.097 


0.47 ±0.17 



Table 1: The predicted abundance of DLA systems, n{z) , for a variety of cosmological models. 
Q,\ = A/(3Hq) where A is the cosmological constant. References for model normalizations 
are listed below. The final row shows observed data and la errors. 



''Bunn & White (1996) 
''White et al. (1996) 

Xen et al. (1994), Kofman, Gnedin, & Bahcall (1993) 
•^Gorksi et al. (1996) 
'^Ma (1996) 

■''Storrie-Lombardi et al. (1996) 




Fig. 2. — Comoving number density of halos with Vc > 37 km s ^, N{z, 37), vs. redshift. 
The models and labehng are as in Figure 0. 
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which contribute to the z = 4 data point have not yet been confirmed by high-resolution 
spectroscopy, however, so the value could decrease with future observations. Hence LCDM 
could be considered acceptable. LCDM2 is an identical model with the lower normalization 
favored by Cen et al. (1994). The lessened power of this model leads to a lower n{z). 
Actually achieving this low normalization while being consistent with COBE would require 
a tilt or a lower h, which in turn would reduce small scale power, further exacerbating the 
problem. 

It is also possible to preserve an infiationary scenario in an open CDM universe. 
OCDM shows a COBE-normalized open model similar to LCDM but with no cosmological 
constant. The open model has a lower mass fiuctuation amplitude and hence a lower n{z) 
than LCDM. An open CDM model with Qq = 0.5, as ~ 0.9 (not shown) skirts the lower 
end of the observational data and might be acceptable within the current uncertainties. 
The amount of small scale power implied by COBE normalization in open models decreases 
rapidly with decreasing (c/. Figure H). For = 0.3 (0CDM2), as is reduced to near 
0.5, and the predicted n{z) is well below the observations. 

The most popular alternative to lowering Qq is to add hot dark matter. This reduces 
as given by the COBE normalization and makes the shape of the power spectrum closer 
to that observed for galaxies. The first incarnations of this mixed dark matter model had 
Qi, = 0.3, which we reproduce here as MDM using the power spectrum and normalization 
given in Ma (1996). This model has strongly suppressed structure formation at high redshift 
(c/. Figure H), and it drastically underproduces DLA absorption. Largely in recognition 
of this problem, = 0.2, shown here as MDM2, has become the favored version. While 
more palatable than MDM, this model still falls 2.8a short for 2 < 2; < 3 and 2.5a below 
observations at z = 4. This problem could possibly be alleviated by further modifications 
such as "anti-tilting" the primeval spectrum to n > 1 or raising the Hubble constant 
to h > 0.5, but these changes would raise as, which is already too high to yield a good 
match to cluster abundances for Qq = 1 (White, Efstathiou, & Frenk 1993; Eke, Cole, & 
Frenk 1996). The most optimistic reading would be that mixed-dark matter can barely 
squeak by the combined constraints of cluster abundances on one hand and DLA's on the 
other. Improvements in the observational data of DLA's and the firming up of theoretical 
predictions by targeted simulations will probably close this window of parameter space. 

These results are consistent with the findings of Mo & Miralda-Escude (1994), in 
which LCDM underproduces absorption at 2 ^ z ^ 3 by roughly a factor of three and 
MDM (with Qi, = 0.3) fails in this range by many orders of magnitude. Kauffmann & 
Chariot (1994) also conclude that standard CDM gives the best match to observations, 
remarking upon the inabihty of MDM models to produce the necessary cold, collapsed 
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gas at higher redshifts. Ma & Bertschinger (1994) find that an Q,, = 0.2 MDM model 
achieves rough agreement with the observations in the range 2.5 < z < 3.5, but they assume 
that the distribution of gas matches the distribution of dark matter in the halos of their 
(dissipationless) simulation. The KWHM simulations show that only a fraction of the gas 
in a virialized halo cools and becomes neutral and that it is distributed primarily in one or 
more small, dense knots; both effects can reduce the damped Lya absorption cross-section 
of a given halo. Ma et al. (1997) have undertaken a new study of DLA absorption in the 
mixed dark matter model, calibrating the Ma & Bertschinger (1994) dissipationless method 
against the KWHM simulation of SCDM, and they also find that the Q^, — 0.2 model 
underproduces the observed DLA absorption. 

The global cosmological parameters in the CCDM and MDM models are the same as 
those in the SCDM model simulated by KWHM, and we are therefore fairly confident that 
our method can accurately cahbrate their predictions of n{z). The TCDM, LCDM, and 
OCDM models, however, adopt different values of ^b, Qq, h, and A, altering the relations 
between overdensity and gas cooling rate and between redshift and the age of the universe. 
Our assumption that a{vc, z) is close to that of the SCDM model is therefore more suspect. 
CKHW find that the DLA cross-section at fixed decreases with time as the gas condenses 
into tighter configurations. If we assumed a similar trend across changes of cosmological 
parameters, then our predictions of n{z) for TCDM, LCDM, and OCDM would go down, 
so in this sense our conclusions about these models are likely to be conservative. However, 
more detailed numerical studies will be required before we can build a definitive case against 
them. 

As a concluding comment, we can generalize a point that we have already made 
regarding the mixed dark matter models. The model that produces the best match to the 
high-redshift DLA data, CCDM, is grossly inconsistent with the mass function of galaxy 
clusters today. The SCDM model, which fits the DLA data acceptably, is at best barely 
compatible with the cluster data, which imply cig ~ 0.55 for = 1 rather than = 0.7 
(White et al. 1993; Eke et al. 1996). Our LCDM model skirts the low end of the DLA 
data and the high end of the cluster data (see Cole et al. 1997). Clearly these are specific 
instances of a general problem, within the broad class of Gaussian/inflationary models 
considered here, since most parameter changes that reduce cluster masses also reduce DLA 
abundances. Between the Scylla of the damped Lja systems and the Charybdis of the 
cluster mass function lies a narrow channel that few cosmological models will pass. 
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